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The perturbative QCD approach, based on Modified Leading Log Approximation 
and Local Parton Hadron Duality, is shown to describe inclusive features of mul- 
tiparticle production in the soft region. Further predictions of this approach are 
discussed. 



1 Introduction 

The study of soft particle production is of particular interest since new infor- 
mation on the physics of confinement can be obtained. An important open 
question in this subject is whether the soft momentum region is smoothly 
connected to the hard region, where the purely perturbative description is 
appropriate, or dramatic changes are present due to the dominance of the 
hadronization process. In this respect, it is worth noticing that the charged 
particle inclusive energy spectra measured in e~^e~ annihilation are surpris- 
ingly close over the whole momentum range (down to small momenta of a few 
hundred MeV) to the predictions of the perturbative QCD approach, based on 
the Modified Leading Log Approximation (MLLAJiln and the notion of Local 
Parton Hadron DuaUty (LPHDpi (see alsol). 

In this paper I discuss: a) how to test the validity of the perturbative 
description in the soft region and b) how to study the sensitivity of experimental 
data to different aspects of this approach, in particular to coherence effects (see 
H for a discussion of the sensitivity of soft particle production to the running of 
the QCD coupling). I focus in particular on the study of the invariant density 
Edn/cfip = dn/dy(PpT in the limit of vanishing rapidity y and transverse 
momentum px or, equivalently, for vanishing momentum \p\ = p, i.e.. 



/o = lim E-^ = - lim E-^ (1) 

y^a,PT^o d-^p 2 p^o d-^p 

where the factor 1/2 takes into account that both hemisphere are added in the 
limit p —> 0. 

The invariant hadronic density in e~^e~ annihilation is shown to approach 
a cTTis-energy independent value at low particles' momentum, in agreement 
with the predictions of the perturbative approach with the inclusion of coher- 
ence. This result would then suggest to extend the validity of the theoretical 



picture down to the very soft region, thus lending new support to the picture of 
LPHD, and it would confirm the relevance of coherence effects in soft particle 
production. To establish this interpretation more definitely, further tests of the 
perturbative picture are proposed. In particular, new predictions for the yield 
of the soft radiation in 3-jet events in e^e" annihilation are presented. A test 
of universality of soft particle production in different reactions corresponding 
to different partonic colour emitters is also discussed. 

The results contained in this paper have been obtained in collaboration 
with Valert A. Khoze and Wolfgang Ochail and have been partially presented 
elsewherdslH. 

2 The invariant density E-^ in e^e~ annihilation 

From the theoretical point of view, the coherence of the gluon radiation requires 
that very soft partons (with long wave length) are emitted by the total colour 
current, totally independently of the internal structure of the jet. According to 
the hypothesis of LPHD, one would then expect that the hadron spectrum at 
low momenturn, v should be nearly independent of the jet energy Sjet, i-e., of 

laJ^sis of the inclusive spectra in t 



irn, V sh( 



he variable logp, 
suggests indeed 



the cms energjciJ'liil. The analysis of the inclusive spectra 
as measured by the TASSOLIj and OPAL Collaboration 
that the soft tail of the spectrum is approximately independent of the cms 
energy in agreement with LPHD expectations Jji order to characterize better 
this phenomenon, it has recently been proposediil to study the invariant density 
Edn/(Pp at very low particle energies of the order of few hundred MeV. The 
important prediction of the perturbative QCD approach is that the invariant 
density approaches in this limit a value independent of cm,s energy. In the 
following, experimental analyses and theoretical predictions for the invariant 
density in the soft region are presented in more details. 

2. 1 Experimental results for low momentum particles 

Let us first consider the invariant density Edn/<fip for all charged particles. 
The angular dependence is actually integrated out in the following, so let 
us call Edn/dPp the quantity j^^- Since no direct measurements of this 
observable have been performed, experimental information can be extracted 
from the published data on the inclusive momentum spectrum, dn/d^p vs. S.p, 
where ^p = log(l/a;p) = log(-y/s/2p), using the relation 

dn E dn 

E—- = - — -—- 2 

d^p Anp^ d£_p 
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Figure 1: a) Invariant density Edn/d?p of charged particles in e^e~ annihilation as a 
function of the particle energy E = ^/p^ + inrvi at m^ = 270 MeV at various cms energies; 
b) the same as in a), but at raii = 138 MeV. 



where E^ — p^ + m\ and rah is an effective particle mass. 

There is not a unique prescription for the choice of this effective mass; on 
one hand, the value of ruh = 270 MeV has been used in t3 to give a good 
description of the moments of the energy spectra with the MLLA formulae. 
On the other hand, since pions dominate in the soft region, the pion mass 
would be a natural value for this parameter. 

Fig. ffl) shows the charged particle invariant density in e+e^ annihilation 
as a function of the particle energy E at different cms energies ranging from 3 
GeV up to LEP-1.5 cms energy (133 GeV^EJ'EJ; two different values of mh 
are used. It is remarkable that the data from all cms energies tend to converge 
in the soft limit. The choice of rrih does not modify the gross features of the 
invariant density, although two different values of the soft limit /g are observed. 
More precisely, let us point out that the LEP data seem to tend to a limiting 
value larger by about 20% than the data at lower cms energies. This may be 
due to the overall systematic effect in the relative normalization of the different 
experiments (in this respect, let us stress that the invariant density is to a 
better approximation energy independent within sets of data collected at the 
same detector) . Alternatively, a possible physical source of energy dependence 
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Figure 2: a) Invariant density Edn/cPp of charged pions in e+e annihilation as a function 
of the particlfi-energy E at various cms energies; b) the same as in a) , but for charged kaons 



and protons I 



The kaon distribution is multipHed by a factor 10 for the sake of clarity. 



is given by particles produced via weak decays. They should indeed be added 
incoherently to the particles produced from the primary quarks and thereby 
could yield a rise of the soft particle spectrum with increasing energy. 

Let us now turn to the invariant density for charged pions, charged kaons 
and protons. Fig. (g) shows the invariant cross section Edn/cPp for tt's, K^s 
and p's as a function of the particle energy E, as extracted from the inc his iv e 
momentum spectra measured at cms energies from 1.6 GeV to 91 Ge\r fr^ 
according to eq. (0). In this case, the simplest choice for the mass scale rrii is 
given by the particle mass itself. In all cases the data tend towards an energy- 
independent limit for E ^ rrii {p ^ 0). Notice that the value of Iq is roughly 
proportional to m^ . 



2.2 Theoretical predictions on the spectrum shape 

The description of the soft region of inclusive energy spectra requires some 
care in the treatment of kinematics. 

Experimental hadronic spectra are usually presented as function of mo- 
mentum p or ^p = log(l/a::p), and they do not have any kinematical boundary 
from above. On the other hand, in the MLLA the partons are treated as mass- 



less but the transverse momentuni is required to be larger than the infrared 
cutoff Qo, such that E — p > pr > Qo- The theoretical predictions are then 
limited from above, as ^ < F = logy/s/2Qo. 

As far as E c^ p ^ Qq, the different kinematical boundaries can be ne- 
glected and the analytical prediction using LPHD is uniquely defined. How- 
ever, the deeper one goes into the soft region, the more sensitive one becomes 
to the different kinematical thresholds. In the soft region, there is not a unique 
procedure for deriving analytical predictions applying LPHD, but one has to 
introduce some additional assumptions and define practical prescriptions. As 
a general rule, one may require that the invariant density Edn/cfip obtained 
at hadron level by one of these prescriptions approaches a constant limit for 
p — *■ as observed experimentally. In the following two possible alternative 
Ansatze are discussed. 

It is important to stress that the most important prediction of the pertur- 
bative QCD approach, i.e., the energy independence of the invariant density 
for particles of low momentum, does not depend on the particular Ansatz one 
uses. This feature will indeed be recovered in both Ansatze we are about to 
discuss. 

Ansatz 1 

IS relate hadron and parton spectra in a single A-jet in the following 



Letus 
wa30ll 



e/-^ = K,E/-f^ . K,D^A(E,Qo. A) (3) 

dph dpp 



with Eh = y^p'l + Qo = Ep = E > Qo, ^E = ^ = \ogQ/E and Kh a normal- 
ization parameter. If hadrons from both hemispheres are added, Kh should be 
replaced by 2Kh. 

Accordingly, the invariant density Edn/d'^p to be compared with the data 
in e'^e~ annihilation is given by: 

dn Kh 2Cp , „ . , 

with E ~ y/s/2e~^ and p^ — E^ — Qq and D{£_, Qo, A) is the inclusive spectrum 
predicted for a single gluon jet. Cp and Nc are the colour factors for quark 
and gluon jets, Cq = Cp and Cg = Nc; the ratio Gp /Nc is needed to obtain 
the spectrum in a single quark-jet and the factor 2 to take both hemispheres 
into account. 

Let us notice that, as far as the theoretical spectrum D{^,Qo,K) linearly 
goes to for p ^- 0, as in MLLA (see below), the invariant density for hadrons 



approaches indeed a finite limit as in ffl) 

with /?2 = ANc/h, b = (lliVc - 2n/)/3, and A = logQo/A, where A is the 
QCD-scale; Nc and nj are the number of colours and of flavours respectively. 

In addition, this prescription is particularly suitable for phenomenological 
analyses of the moments of the inclusive energy spectruniij. In this respect, 
let us notice that a good phenomenological description of the charged particle 
inclusive spectrum has been obtained with Qo — 270 MeV; this value of Qq 
will then be adopted in the theoretical calculations. 

To obtain the theoretical prediction for the invariant density Edn/dPp 
within this scheme, let us then obtain an explicit solution for the spectrum 
£'(^, QojA) in the soft limit. Consider first the Double Log Approximation 
(DLA), in which energy conservation is neglected and only the leading sin- 
gularities in the parton splitting functions are kept. The evolution equation 
of the inclusive energy distribution of partons p originating from a primary 
parton A is given bytll: 

D\{tY) = 5\5{0 + j^ dej^ dy'^lliy' + ODHe.v' + O (6) 

where S, — log(l/2;) = \og{Q/E) and Y ~ ^og{Q/Qo) with E the particle 
energy and Q the jet virtuality {Q — PO for a jet of primary momentum P 
and half opening angle 0); Ca is the respective colour factor for quark and 
gluon jets; 70 denotes the anomalous dimension of multiplicity and is related 
to the QCD running coupling by Jq = iNcasf^n = (3^ / \og{pT / ^) ■ 
Eq. (0) can be solved iteratively; with two iterations one gets: 

D^Ai^,Y) = S^J{0 + ^P'log^l + ^^ (7) 

^ (3^J^~^dTlogil + l)log{l + ^j' 



log(l + ^) 

The term of order /3^, corresponding to a single gluon emission, yields the 
leading contribution for E — > Qq. It is worth noting that this term does not 
depend on the cms energy and it is proportional to the colour charge factor 
of the primary parton. Therefore, we have explicitly found the features of 
soft radiation which we have qualitatively expected. In addition, the second 

6 



E— 
iPp 

[GeV-^ 



0.1 




MLLA (A = 0.01) 
DLA (A = 0.01) 
Limiting Spectrum (A : 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0. 



E[GeV] 



Figure 3: Invariant density Edn/ifip as a function of the particle energy E for Qq = 270 
MeV. Predictions at cm^ energies of ^/s = 3 GeV (lower three curves) and 91 QeV (upper 
three curves) using eq. (kh with Dg computed in MLLA (eqs. (pl,p|)), DLA (eq. (pi)) and the 
Limiting Spectrum (normalization Kf^ = 0.45, 0.45 and 1.125 respectively). 



iteration provides us with a energy-dependent term, which allows to describe 
the rise of the spectrum with increasing ^/s at particles' energies around 1 
GeV. Let us also point out that the ^-spectrum vanishes in the soft limit 

^^Y {E^Qo) as 



D^i^^Y) 



Y-^r^logE/Qo^E-^ 



(8) 



In order to get an improved description at nonasymptotic cms energies, it 
is important to use the MLLAna, which takes into account the exact form of 
the parton splitting functions in the evolution equation. The MLLA solution 
is in general much more involved than the DLA one, but in the soft limit the 
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Figure 4: Invariant density Edn/d^p of charged particles in e+e annihilation as a function 

of the particle energy E = x/p^ + Qq at Qo = 270 MeV. Data points at various cms energies 

from SLAG, TASSO and TOPAZ Collaborations, LEP-1 and LEP- 1.00113 are compared 
to MLLA predictions (A = 0.01, Ky, = 0.45). 



expression simplifies and one gets a simple exponential suppression factor: 



D{^,Y, X)\mlla ^ D{^,Y, X)\DLAexp 



a / ^^{y)/{mc)dy 



(9) 



with a = llNc/S + 2nf/3N^. This solution is valid in the soft limit only, when 
the difference [7o(C)~7o(^)] i^ small compared to j^D. As a consistency check 
of eq. (P) , the MLLA with fixed coupling has been explicitly solved and eq. (M 
has been found to be exactly satisfied in the full energy region. 

In order to estimate to what extent the iterative solution gives a good 
description of the inclusive spectrum, one can consider the inclusive spectrum 
D{(^,Qo,A) in the limiting case Qq — A. In this case, the solution, called 
Limiting Spectrum, is explicitly knowra and one can avoid the iterative pro- 



cedure. Notice that the Limiting Speetrum does not have the same behaviour 
near the boundary ^ ^ y as the full solution which follows from the inte- 
gral equation (g) with eq. (g), as the Limiting Spectrum goes to a constant 
as -B — ^ Qq. The invariant density computed from the Limiting Spectrum 
according to eq. (||) shows then a singularity in the soft limit, p ^ 0. 

To illustrate the above analytical results, the predictions for the invariant 
density Edn/d^p in the low particle energy region at the two cms energies 
of 3 and 91 GeV, obtained via eq. (|) in DLA (eq. @)) and MLLA (eq. (§) 
with A = 0.01 and in the Limiting Spectrum case are shown in Fig. (|^). The 
normalization of the limiting spectrum is as in fits to e+e^ annihilation; the 
normalizations of the DLA and MLLA curves are chosen to approach the lim- 
iting spectrum for energies E > 0.5 GeV. 

Both in DLA and in MLLA the invariant density approaches an energy 
independent value in the soft limit ^ ^ Y {E ^ Qq). This originates from 
the soft gluon emission contribution of order a^ which is determined by the 
total colour charge of the primary partons due to the colour coherence. In 
this limit the MLLA converges towards the DLA. The MLLA solution is in 
good agreement with the Limiting Spectrum solution, thus showing that the 
iterative solution includes indeed the dominant contribution in the soft region. 

A quantitative comparison of the theoretical predictions of MLLA with 
experimental data for the invariant density of all charged particles is shown in 
Fig. (Q). A good agreement is visible, both in the energy independence of the 
soft limit Iq and in the energy dependence of the invariant density at particle 
energies above 500 MeV. 

It is also possible to compute the prediclions of MLLA with fixed coupling. 
In this case a full solution can be obtainedi^a. It is interesting to notice that 
the model with fixed coupling cannot describe the steep slope at very low 
momentum, thus suggesting that the running of the coupling is relevant in this 
region. Further studies of models, where the coupling is frozen below a certain 
momentum threshold, are in progress (see alsoli^fl for further discussions of 
the effect of the running of the coupling in inclusive energy spectra) . 

Ansatz 2 

Let us consider the alternative identificationiMI: 

|i^(|)'A,„.(&,y) (10) 

where I?iim(C-Ej ^) indicates the Limiting Spectrum solution. With this Ansatz, 

9 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0. 



E[GeV] 



Figure 5: Invariant density Edn/d^p as a function of the particle energy E = Wp'^ + Qq 
with Qo = 138 MeV predicted at various cms energies using eq. (0) {Kfi = 1.125). 



the invariant density for hadrons in e+e annihilation is given by: 



dn Kn 2Cp 



(11) 



The preferred value for Qo = A is in this case the pion mass, as natural in the 
soft region, where most produced particles are indeed pions. 

Let us notice that as p ^ 0, the particle energy E goes to the finite value 
Qo and Dum goes to a constant, i^fij^. One then recovers again a finit value 
for the soft limit 



lo 



Kh 2Cf 
4^gg Nc '™ 



(12) 



The invariant densities predicted with this method at different cms energies 
are shown in Fig. (||). It is important to stress that also with this second 
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Figure 6: a): invariant density Edn/d p of charged particles in e+e~ annihilation as a 
function of the particle energy E = Wp^ + m| with m^ = 138 MeV at various cms energies 
compared with the predictions of eq. (hlh with Qo = 138 MeV and Kf^ = 1.125; b): same 
as in a), but for charged pions. 



Ansatz all curves tend to converge to a common value independent of cms 
energy at very low particles' energies. 

This approach has been applied both to the invariant density of all charged 
particles and to the invariant density of pions. Its predictions at different cms 
energies are compared in Fig. (g) with experimental data for charged particles 
and charged pions. At large cms energies, a good agreement is found for 
charged pions in the whole particle energy region and for charged particles 
in the low particle energy region. However, this prescription cannot properly 
describe data at low cms energies, i.e., it does not reproduce correctly the cms 
energy dependence of the data. 

3 Further tests of LPHD in 3-jet events in e+e annihilation. 

In view of the success of the LPHD picture for describing the invariant density 
of charged particles in e~^e~ annihilation, it is interesting to look for other 
independent tests of this picture. In this section, a new analysis which could 
shed new light on the validity of the perturbative description of the soft region 
is proposed. 

Consider 3-jet events, measure the invariant density of particles produced 
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into a cone perpendicular to the production plane for different angular con- 
figurations of the three jets and extract then the limiting value /q for each 
configuration. To get an absolute normalization, look also at the radiation 
into a cone of the same angle again perpendicular to the primary qq direction, 
but in 2-jet qq events and build the ratio 

The proposed observable seems at a first sight similar to the study of the 
string phenomenon, which concerns the radiation in the plane of 3-jet events; 
the advantage of studying the radiation in the transverse direction lies in the 
fact that in this case one can avoid the integration over the A;t > Qo boundaries 
along the jets. Moreover, the new observable refers to soft production only, and 
not to the total particle flow integrated over momentum, as in the discussion 
of the string effect. For the particle flow, the original angular flow should be 
convoluted with a cascading factoHil, which, however, cancels out when one 
looks at ratios. The same predictions for the ratio i?_L should then in principle 
be valid both for the ratio of the soft particle yields, /q, and for the ratio of 
the multiplicity flows. However, it would be useful to perform independently 
both measurements. 

Let us now consider the theoretical prediction for the above observable 
in the framework of the perturbative approach. The calculation only includes 
the contribution of soft gluon bremsstrahlung emission up to order ag, since, 
as seen in the study of the invariant density in e+e^ annihilation, this term 
should dominate in the soft limit. 

The formula for the soft radiation into arbitrary direction n from a qq 
antennapointing in directions n^ and Uj is given by the usual bremsstrahlung 
formul£iJ: 



dp ,„ as 
—dVL^—— 



dN,^^^dnn—^W<i^{n) , M/«nn) = 2Cf(zj) (14) 



with {ij) ~ aij / {aiaj), Uij = (1 — fiiTij) and Oj = (1 — nfii). For the radiation 
perpendicular to the primary partons (ij) — aij = 1 — cos Oij, with Oij the 
relative angle between the primary partons i and j. 

The soft gluon radiation in a 3-jet event is given as in eq. (14h with the 
angular factor 

1 



W^^^ {ft) = Nc m) + ifq) - -^ m] (15) 
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R± 


i?i (large Nc) 


^99 =^ - Q99 


1 


1 


(collinear or soft gluons) 






99 ^ 99 ^ ^^ 


1.59 


1.5 


(Mercedes) 






6gg = Qgq — TT 


ff = 2.25 


2 


(gg antiparallel to g) 







Table 1: Prediction Q) for the ratio iJx = dN"^^ /dNf and the large-Afc-limit (|l8|) for 
different angular configurations of the qqg events (Sqq = 2it — &gq — &gq). 



The soft gluon radiation for a qq pair with relative angle Qqq is given again 
as in eq. (14), but with the angular factor 



WfiQqq) = 2Cf{1 - cose,,-) . 



(16) 



which, for 2-jet events in the rest frame, gives W]j'{7t) — ACp- 

Correspondingly, the ratio R± of the soft particle yield in 3-jet events to 
that of 2-jet events in their own rest frame is given by 



Ra 



dm 



dN 



99 



^^[2 -cose,. 



cose 



99- ^^(1- cose,,-)] (17) 



m 



Theoretical predictions for R± are presented in Table (|l|) for three different an- 
gular configurations; notice that in the two extreme cases, the soft or collinear 
primary gluon emission and the parallel qq configuration, the expected lim- 
iting values are correctly recovered. The particularly simple and interesting 
situation of Mercedes-type events, where no jet identification is necessary for 
the above measurement, is also shown. The table also shows the results ob- 
tained in the large- A'^c approximation, in_azliich the qqg event is treated as a 
superposition of two qq dipoles (see, e.g.,l££l). In this case the expression for 
the ratio R± simply becomes: 



Ra 



1, 



[2 — cos 6(1 



cos 60 



(18) 



The differe: 
investigatei 



among the full theory and the large-A^;;; limit can also be 
121 by studying the production rate in 3-jet events normalized to 
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the sum of rates from the corresponding 2-jet events (dipoles) with opening 
angle Qgq and Qgq respectively, rates which can be extracted experimentally 
from the analysis of qq^ events. 

4 Test of universality in different reactions 

Let us look at the soft limit of the invariant density 



/o = lim E-^ = - lim E-^ (19) 

in different reactions. This study provides a very direct test of universality in 
soft particle production; one could indeed directly measure whether the soft 
particle production is universal, i.e., it is a purely hadronic quantity not related 
to the underlying partonic processes, or the intensity /q depends on the colour 
topology of the primary active partons in the collisions process, as predicted 
in our perturbative approach. For instance, in case of quark exchange the two 
outgoing jets originate from colour triplet charges and Iq should be as in e+e~ 
annihilation; on the contrary, in case of giuon exchange Iq should be about 
twice as large (Nc/Cp)- If data would show these features, a further direct 
support of the dual description of soft particle production in terms of the QCD 
bremsstrahlung would be achieved. 

Unfortunately, an absolute prediction of Iq is not possible, since it depends 
on the normalization factor and the cut-off parameter Qq. One should then 
exploit the energy independence of the soft limit Iq in e+e^ annihilation and 
use this value to set a standard scale to be used for the comparison with other 
processes. 

4-1 Quark exchange processes 

In case the process goes through quark exchange, i.e., via colour triplet ex- 
change, the soft production intensity Iq should be the same as in e^e^ anni- 
hilation. 

A possible realization of quark exchange process is the process 77 —> qq —^ 
2-jets with either the virtuality Q^ of t he initial photon or the scattering angle 
photon-jet sufficiently large (see, e.g. c3). Another example is deep inelastic 
scattering at large Q^. In this case the current fragments in the Breit frame 
are expected to have the same characteristics as the qua rk f ragments in one 
hemisphere of e+e^ annihilation (for a QCD analysis, se£3), and this is in- 
deed observed for not too small Q^ c2l. A third case is the production of two 
W's which decay in the fully hadronic channel. In this case, since each of 
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the two W's decays into 2 jets, the limiting soft particle yield should he twice 
the yield in e+e~ — > qq , provided the interconnection phenomenaCd are ne- 
glected. Alternatively, deviations from the factor 2 could give an estimate of 
the importance of these effects. 

4-2 Gluon exchange processes 

In case the process goes through gluon exchange, i.e., via colour octet exchange, 
the soft production intensity Jq should be roughly twice as large {Nc/Cp) as 
in e+e^ annihilation. The same factor appears in the theoretical predictions 
for the ratio of the average multiplicity in gluon vs. quark jetcll; experimental 
analyses of this ratio have obtained values considerably below this limit (see 
e.g. u). Since in our case we are not integrating over momentum but we 
are simply looking at the soft radiation, where some effects, like for instance 
energy-momentum constraints, are expected to be small, the effects of the 
different colour charges could be more pronounced than in the study of the 
multiplicity ratio. 

The realization of a gg colour singlet final state is not easy experimentally. 
It may become available at future colliders at higher energies through the 
process 77 — > gg. An approximate realization of a colour octet antenna has 
been recently obtained in e~^e~ —^ qqg, by selecJing a configuration with the 
gluon recoiling against a quasi-coUinear qq paii» (for a recent study of such 
events, seeLiJ). A process mediated by colour octet exchange is expected to 
occur also in DIS for Q^ > few GeV^ and small Bjorken x, via the diagrams 
of photon gluon fusion. Another example is hadron-hadron collision with a 
particle or jet at moderate pr (> 1-2 GeV) at small angles so that the overall 
2-jet structure is maintained. 

4-3 Soft collisions (minimum bias events) 

These processes (with initial hadrons or real photons) are not so well under- 
stood theoretically as the hard ones but it rnight be plausible to extrapolate 
the gluon exchange process towards small p^eil. Accordingly, the soft radia- 
tion is expected to be twice as large as in e^e^ annihilation. For these reac- 
tions, some experimental information is already available. Data at ISR energies 
^t the soft limit Iq is similar to e+e" annihilation at a similar effective 
On the other hand, Iq roughly doubles when going from y^ '^ 20 
to Y^ = 900 GeV at the collidet^. If additional incoherent sources 
can be excluded, such behaviour could indicate the growing importance of one- 
gluon exchange expected from the perturbative picture. Then a saturation at 
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Iq^/Iq ^ ~ 2 and no additional increase for a semihard pt trigger would be 
expected. However, the rise of /g '' ^-t collider enereies could also result from 
the incoherent multiple collisions of partons (e.g. 123) which has recently been 
postulated in 7p collision£ll. This unclear situation certainly deserves further 
studies. 

^.J^ Rapidity dependence of Iq in DIS 

According to the previous discussion, different kinematical conditions in jp 
collisions select different elementary subprocesses and should then give rise to 
different yields /q's. For example, by decreasing Q at fixed hadronic energy 
W one goes from quark exchange process to a gluon exchange process and 
finally to a soft process. We have seen that different reference frames, for 
instance the Breit frame, can be useful to select a particular situation. An 
alternative possibility is provided by the study of the soft particle yield /q at 
different values of rapidity in one reference frame, for instance the cms frame. 
Let us consider then the yield Io{y) — ^j-jr — | _,g as a function of rapidity 
y, measured in the cms frame. One would then expect the existence of a 
"quark plateau" of Io{y) in the current region of length Ay « \og{Q/m) and 
height equal to e^e^ annihilation near the rapidity corresponding to the Breit 
frame, i.e., in the direction of the incoming photon, y Breit', moving towards 
the central region, a transition to a "gluon plateau" in the complementary 
region is expected to occur. If sufficiently hard gluons are exchanged in the 
process, Iq should then increase but still remain smaller than twice the value 
in e~^e~ annihilation. Approaching the fragmentation region, a transition to 
soft interactions should be observed and a decreasing of Iq is likely to occur. 
In general, the soft particle yield Iq as a function of the rapidity should show 
a step like behaviour, pointing out the regions where the underlying process 
is dominated by quark or gluon exchange; a direct indicator of the underlying 
mechanism would be provided directly by the height of the plateau. 

5 Conclusions 

The analytical perturbative approach to multiparticle production, based on the 
Modified Leading Log Approximation and Local Parton Hadron Duality has 
been shown to successfully describe data on inclusive energy spectra of charged 
particles in QCD jets. In order to investigate the limitations of this picture 
and to point out non perturbative effects, this approach has been applied in 
the soft region. 

The invariant densities Edn/d^p of all charged particles and identified 
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particles in e+e^ annihilation arc found to be approximately energy indepen- 
dent at low particle energy over a very large range of cms energies. The same 
behaviour is expected in the perturbative approach, as a consequence of the co- 
herence of the soft gluon radiation from all emitters. With proper additional 
assumptions for the treatment of mass effects, the perturbative predictions 
reproduce the data also quantitatively, thus extending the phenomenological 
success of the perturbative approach down to the soft region. This result sug- 
gests that the production of hadrons in the soft region, which is known to pro- 
ceed through many resonance channels, can be simply parametrized through 
a parton cascade pushed down to a small scale of a few 100 MeV. 

The study of transverse radiation in multijet events in e+e^ annihilation 
has been proposed to test the sensitivity of the soft particle production to the 
effective colour charge of the primary emitters. Predictions of the perturbative 
approach for the soft particle production in different reactions have also been 
presented. The perturbative approach predicts a breakdown of universality 
among different reactions, since the soft particle yield reflects in this picture 
the properties of the underlying partonic process. 
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